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SUMMARY 

A method is presented for the determination of the contour of 
disks, typified by those of aircraft gas turbines, to incorporate 
arbitrary elastic-stress distributions resulting from either cen- 
trifugal or combined centrifugal and thermal effects. The specified 
stresses may be radial, tangential, or any combination of the two. 
For example, the octahedral shear stress may be specified at each 
radius inasmuch as this stress Is generally used as the criterion 
of the occurrence of plastic flow. Use is made of the finite- 
difference approach in sol7lng the stress equations, the amount of 
computation necessary in the evolution of a design being greatly 
reduced by the Judicious selection of point stations and by the aid 
of a design chart. Use of the charts and of a preselected schedule 
of point stations is also applied to the direct problem of finding 
the elastic- and plastic-stress distribution in disks of a given 
design, thereby effecting a. great reduction in the amount of cal- 
culation compared with previously published methods . Illustrative 
examples are presented to show computational procedures in the 
determination of a new design and in analyzing an existing design 
for elastic stress and for stresses resulting from plastic flow. 


INTRODUCTION 

The important role of the disk in the gas turbine has directed 
much attention to the analysis of stresses in rotating disks with 
temperature gradient. Early in the development of stress analysis 
of steam turbines, Stodola (reference l) presented the basic equa- 
tion for stresses in such disks, as well as a limited approach 
towards the solution. More recently, in connection with the develop- 
ment of the gas turbine for aircraft propulsion, the problem of 
determining thermal stresses resulting from temperature gradient 
has been actively investigated. Several methods have been developed 
(references 2 to 5) that can be used to calculate the stresses in a 
given disk due to differential expansion and to variation in physi- 
cal properties of various parts of the disk operating at different 
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temperatures. Thus, for disks already designed, it is readily pos- 
sible to determine with a high degree of accuracy the centrifugal 
and thermal stresses resulting from a given speed and temperature 
distribution. 

/ 

Little attention has, however, been directed towards the prob- 
lem of designing the profile of a disk, which, when operated at a 
given speed and temperature distribution, will incorporate a given 
total-stress system - centrifugal plus thermal. A direct method of 
design for prescribed total-stress systems was developed at the 
NACA Lewis laboratory. Although the design method can readily be 
extended to include plastic flow, in the interest of simplicity 
only elastic stresses are considered. 

The specified stresses may be either radial or tangential 
values, as is common in engineering practice, or the specification 
may refer to any combination of these stresses. For example, spec- 
ification of the octahedral shear stress at each radius may be 
desired Inasmuch as this stress is frequently used as a criterion 
of the occurrence of plastic flow. The method is based on the 
finite-difference equations derived in reference 3 and makes use 
of several simple design charts, presented herein, that facilitate 
the calculations. An illustration is presented to trace the actual 
steps that may be followed in evolving a design. 

A simplified method of analyzing an existing disk for cen- 
trifugal and thermal stress is also described. The method is espe- . 
daily useful in checking the stresses in a disk designed for pre- 
scribed stress distributions, particularly if minor deviations from 
the design profile are introduced for practical reasons, or if creep 
or plastic flow, not considered during the design stage, is likely 
to occur. Applications to arbitrary disks are equally, valid . 
Essentially, the method is an adaptation of the finite -difference 
approach of references 3 and 4, the principal modifications lying 
in the use of a preselected distribution of point stations and the 
addition of several alinement charts to facilitate the computations. 

Although, the methods presented are of special advantage, in 
calculations involving thermal stresses, they are suitable for disks 
subjected solely to centrifugal effects. Simple equations may be 
obtained both for design or for stress analysis of rotating disks 
by setting ail terms involving temperature equal to zero. The 
subject of design for prescribed centrifugal-stress distribution 
has likewise received only limited attention (examples of which 
may be found in references 6 and 7); the present method may there- 
fore satisfy a need in this field of design. 
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DIHECT METHOD OF DISK DESIGN 


Equations. - The equations for stress analysis of rotating 
disks are (reference 3) 

TT (rho^-hOfc+p^hr 2 = 0 (l) 


and 


d_£t\ 

dr \E / 


_d_ 

dr 



+ 



(i-Hi)(a r -a t ) 

Er 


= 0 


(All symbols are defined in the appendix. ) 


( 2 ) 


Equation (2) is a statement of the relation that must exist 
between the radial and tangential stress distributions, regardless 
of the shape of the disk, and therefore does not contain the disk 
thickness h. Once the stresses have been chosen to be compatible 
in accordance with equation (2), h can be determined from 
equation (l). 


Solution of compatibility equation for specified stresses. - 
A finite-difference method for the solution of equations (1) and 
(2) is presented in reference 3. Choosing discrete point stations 
along a radius permits the stresses at the nth station to be 
expressed in terms of the stresses at the (n-l)st station. Equa- 
tion (2) in finite -difference form becomes 


c 'n^n^'n^n = ^'n^n-l^V^n-r^’n 


(3) 


The coefficients C' n , D' n , and so forth are, in general, com- 
plicated expressions involving the radii r n and r n-1 at the 

nth and (n-l)st stations as well as the properties of the disk 
material at the two stations. A great simplification in the expres- 
sions results if the radii r n and r^-^ are not chosen at random, 

but are chosen according to a preselected schedule, for example, 
according to the relation 


r n“ r n-l 


0.200 


If, in addition, Poisson's ratio u is assumed constant and equal 
to 0.333, equation (3) reduces to 
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1.1330^-0.4660^ = ^■(0.833a t#n . 1 -0.166a r ^. 1 )-E a H' n = ^ 


(4) 


Equation (4) is useful in the central region of the disk where 
values of r n are relatively low. In order to obtain reasonable 

spacing of the point stations in the rim region where the radii are 
large, it is desirable to choose the spacing according to the 
relation 


r n' r n-l 


.0.050 


Equation (3) then becomes 


®n 


1.033a t)n -0.366O r>n = ^(0.9650.^-0.2980, ^>-1^ . ^ 


(5) 

As an example, equation (4) is considered. If ^ r>n _i and 
°t n-1 are a ^ re& dy established and if the temperatures at the nth 
and (n-l)st station are known, the value of can be determined 
then the only unknowns are n and a t Q . The equation is linear 
in O r ^ n and O t n so that whatever are the actual values, a plot 
of o r n against n results in a straight line defined by 

equation (4). In order to establish the individual values of 
O r n and a t n , another condition must be specified. In some 

designs, specification of either the radial or the tangential stress 
may be desired; hence, the unspecified stress may then be directly 
determined from equation (4) or from the corresponding straight-line 
plot. A disadvantage of starting a design by specifying the radial 
or tangential stress is that the unspecified stress may turn out too 
high or in scans other way unfavorable. It may be more desirable to 
specify the octahedral shear stress 0 e , which is defined by the 
equation 


'e.n 




n "°r,n a t,n +a t,n 


( 6 ) 


This stress is usually used as the criterion for plastic-flow cal- 
culations (reference 4) and Is therefore a logical design parameter. 



NACA TN 1957 


5 


Furthermore, specification of the octahedral shear stress Insures 
that neither the radial nor the tangential stress will be unex- 
pectedly high. A plot of O r ^ n against for a given value 

of O- _ results In an ellipse Inclined at 45° to the coordinate 

0 y U 

axes. The intersection of the line defined by equation (4) and the 
ellipse of equation (6) defines the specific values of O r Q and 


G + _. Although, In general, a straight line Intersects an ellipse 

W yH 

at two points, the proper point for design use Is usually obvious 
from the physical nature of the problem. 


Design chart. - A design chart based on equations (4) and (6) 
is presented In figure 1(a). The straight lines correspond to dif- 
ferent values of 2^ In equation (4) and the ellipses correspond 

to different values of octahedral shear stress o Q ^ n . This chart 

Is used to calculate the stresses at the nth station after the 
stresses have been determined at the (n-l)st station. The value 
of ^ is calculated by equation (4) from conditions at the 
(n-l)st station and the temperature at the nth station, thereby 
establishing the line along which G r>n and mU8 ‘ fc 

Intersection of this line with the ellipse corresponding to the 
assigned value of octahedral shear stress at the nth station 
establishes C r Q and Q . Thus, starting at a station near the 

center of the dish (or for a dish with a central hole, at the edge 
of the hole) where the radial and tangential stresses bear known 
relations to the octahedral shear stress mahes it possible to pro- 
ceed with the determination of stresses at successive stations of 
Increasing radius . 


The corresponding chart for equation (5) Is figure 1(b), which 
Is used in determining the stress conditions In the rim region, 
whereas figure 1(a) Is used for determining the stress conditions 
In the central region of the dish. 

Determination of dish profile. - When the stresses have been 
determined, the shape of the dish. can be obtained from the finite- 
difference equivalent of equation (l). For random values of r n 

and r n-1 , the equation assumes the form 


Vi 

V 


(^^ r ’ n ' 0t ’ D+0n ^ rn2 

2r n-l\_ _ 2 2 

-P n-1 (,) r n-l 


( 7 ) 
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r — r *i 

For — — Sli s 0.200 and Pq ■ Pn-i = constant density, equa- 


tion (7) becomes 


Vi 

*n 


10.000 -or +P r ^o 2 
r,n t,n n n 

2 2 

8.000 ,+C . -0.640p r CO 

r,n-l t,n-l n n 


(8) 


X* “X* -I 

and for — — 2l± = 0.050 and p Q = P n _i = constant density, eq.ua- 
r n 

tion (7) becomes 


Vl . *0-°0«V,n-°t > n^ r n “ (9) 

^ M-«» r ,n-l +<> t,n-l- 0 - 903 ‘ > n r n 2 “ 2 

Because equations (8) and (9) give only ratios of thickness 
and not specific thicknesses, a starting point for calculating actual 
thicknesses must be obtained from other considerations, for example, 
those at the rim. For disks with the familiar fir-tree blade 
attachments, one approach is as follows: 

With reference to figure 2, if sufficient clearance is provided 
in the fastening between the blades and the disk to prevent circum- 
ferential tightening, the stress in the serrated region of the disk, 
neglecting stress concentrations in the teeth, is essentially uni- 
axial. The radial stress in the disk at the base of the serrations 
(at radius r^) . is the sum of the centrifugal stresses produced by 
the blades and by the serrated regions of the rim. 

The total centrifugal force due to airfoil sections of blades 
is 


HiV < 10 > 

The serrated bases of the blades and the base -retaining sections of 
the disk may be considered as a discontinuous ring of material sub- 
jected to centrifugal effect, the serrations serving only to elim- 
inate tangential stresses. Assuming that at the rim the thickness 
hd (fig- 2) is specified from considerations of blade design and 

that the thickness tapers linearly to a value h^ (as yet to be 

determined) at the base of the serrated section of the disk, then 
approximately : 
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Total centrifugal force on serrated sections of blades and 
disk is 

p 

(speed 6 ) (volume) (average density) (average radius) 


or 


(i) 2 (r.-r,) ( W Zn 


(ID 


The radial stress at the base of the slots is equal to the radial 
stress O r b at the rim of the continuous section of the disk and 

presumably has already been determined from the disk design. Hence, 
Total resisting force at the rim of continuous disk 


= (stress) (area) 
= °r,b 2 * r b h b 


( 12 ) 


When the sum of equations (10) and (11) is equated to equation (12), 


“%iV W( V r t> )2( V r b ,h d 

h b = (13) 

8rtr b a r,b- np ^ (r d +r b )2(r d“ r b } 


When h^ has been determined, the thicknesses at the remaining 
stations in the disk may be obtained by successive multiplication 

by the ratio An application of the method is Illustrated by 

*n 

the following example. 


ILLUSTRATIVE DESIGN 

Requirements . - The assumed problem is to design a disk of 
9-lnch radius to carry blades that require a 1.9-inch base thick- 
ness (h^ =1.9). In order to retain the blades, a serrated region 
of 1-inch radial depth is necessary (r^, = 9-1 = 8 in.). Based on 

measurements of similar disks, the temperature distribution is 
assumed to be as shown in figure 3, and as an initial design, the 
equivalent stress at each location is assumed to be equal to the 
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elastic limit of the material at the operating temperature. The 
disk is to operate at a speed of 11,850 rpm and the density of the 
disk and blade material is 0.29 pound per cubic inch. 


.Establishment of compatible stresses. - The calculations are 
shown in table I. The disk is divided into 18 stations. Column 1 
indicates the radius at each station as a fraction of the radius 
of the continuous section of the disk. These values have been 


chosen to yield 


with 


r n~ r n-l 

r n 


r -r i 

9 stations with =r— = 0.200 and 9 stations 

n 

0.050. The actual value of r n at each station is 


obtained by multiplying the values of column 1 by the disk radius, 
which in this case is 8 inches. The radii are listed in column 2. 


The difference between the actual temperature and the tem- 
perature at which there is no thermal stress is listed in column 3. 
Because in this case the stress-free condition is at a room tem- 
perature of 70° F, the AT n term of column 3 is obtained by sub- 
tracting 70° F from the actual temperatures in figure 3. 

The elastic moduli and the coefficients of thermal expansion 
at the operating temperatures of the various stations are listed in 
columns 4 and 5, respectively. 

The values V T n' H’ n = VV’h-A-P “‘'Vn *** 
then calculated as shown in columns 6, 7, and 8, respectively. 

The quantity r n 2 p n 0^ at each station is then calculated, as 
shown in column 9. Column 10 lists values of the ratio 
needed for subsequent calculations. ®n-l 


Columns 11 and 12 are temporarily bypassed, and in column 13 
there are entered the desired equivalent stresses at each station. 

In this design, the octahedral shear stresses are taken equal to 
the elastic limits of the material at the operating temperatures, 
which for the temperatures of figure 3 are shown in column 13. 

The essential stress determinations are then effected by simul- 
taneous operation on columns 11 to 15. At station a, the radial 
Anfl tangential stresses are both equal to O e and are therefore 

both entered as 76,500 pounds per square inch in columns 14 and 15. 
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The 86 stresses are used to determine the value of Z at station 2 

equal to 49,600, as listed in column 12. Reference to figure 1(a) 

then shove that the line Z = 49,600 Intersects the ellipse 

C = 76,300 at „ * 76,000 and at <J. = 75,000, which are 

e J r,n 7 

then listed in columns 14 and 15. In the same manner, these values 
are used to determine the stresses at station 3 and progressively 
at subsequent stations. For example, at station 6, Z = 32,200 
and O e = 75,000. Intersection of the corresponding straight line 

and ellipse yields, in figure 1(a), C r n = 83,000 and 

O. = 62,000. 

r n~ r n-l 

At station 10, there is a change to — — ■ = 0.050. The only 

r n 

effect is to change the coefficients in column 11 and to Introduce 
the use of figure 1(b) Instead of 1(a). At station 10, for example, 
Z > 13,800 and o e a 73,000. Intersection of the corresponding 

straight line and ellipse in figure 1(b) yields = 83,500 

and a t n a 42,500, which are listed in columns 14 and 15, and are 

subsequently used to obtain the value for Z at station 11. Thus, 
by successive and simultaneous operation on the various columns, 
values of <3 r>n and may be obtained for each station that 

are consistent with the compatibility equation and combine to yield 
an octahedral shear stress at each station of any desired value. 

(If the Imposed requirements are impossible to attain, the necessary 
straight line and ellipse will not intersect. Thus, for example, 
were it arbitrarily required to achieve an octahedral shear stress 
at the rim of 40,000 pounds per square inch while maintaining the 
equivalent stresses at all other stations at the values listed in 
column 13, the line Z = 60,900 would not Intersect the ellipse 
O e = 40,000. The indication is merely that no physical disk can 

satisfy the requirements and that the requirements must be changed.) 

Determination of disk profile. - Having established compatible 
stresses makes it possible to obtain thickness ratios in accordance 
with equations (7) and (8). Before proceeding, however, examination 
of the radial stress at the rim 0 r ^ should be made in order tb 

determine whether or not it will yield a reasonable value of thick- 
ness hjj in accordance with equation (13). Substitution of the 


values 
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N 

= 54 

= 27,000 

A bl 

= 0.625 

P 0 ^ = 1150 

r b 

= 8.00 

r d = 9.00 


h d = 1.9 

In equation (13) results In the relation hetvean hf, and O^t, 
shown in figure 4. Obviously, from this figure, values of a r ^ 

below 18,000 pounds per square inch vlll yield inordinately large 
thickness values of h b , whereas values of C r ^ above 

40,000 pounds per square Inch will yield Inordinately low values 
while overstressing the rim material, which is operating at a high 
temperature. The allowable radial rim stress can be ascertained 
from a detailed analysis of creep strength and stress-rupture values 
of the rim material for the desired life of the wheel, but in any 
case, figure 4 indicates that the reasonable range for is 

between 18,000 and 40,000 pounds per square inch. 

Obviously, the value of <J r b in table I obtained in the first 

determination ( 0 _ >. = approximately 0 ) does not lie within the 

r, d 

aforementioned reasonable limits. The reason for the very low value 
of a r b , however, is the designation of the low value of 

O , =58,000 pounds per square inch at rim of continuous section 
of disk. The compatibility condition demands only that Cf^.^ and 
q . lie along the line Z = 60,900 (column 12, table I). By 
modification of the demand on O el) , the individual values of 

and O r ^ can be moved to other locations along Z = 60,900. 

Thus, for example, by a change in the requirements on <7 e ^ to 
63,500, 64,000, 67,700, and 72,900 pounds per square inch, the values 
of a r 15 can be successively moved to 18,000, 20,000, 30,000, and 

40,000 pounds per square inch with corresponding changes in O^b’ 

all these values are given in table I. Physically, the significance 
of this result is the impracticability of attempting to find a disk 
for which stresses at the rim will not exceed the elastic limit if 
all the other regions are to operate exactly at the elastic limit. 
Allowing the rim to operate somewhat above the elastic limit pezmlts 
the remainder of the disk to be kept at the elastic limit. Under 
unusual conditions, modifying several stations near the rim may be 
necessary. 
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The values of 


Vl 


are now computed according to equations (8) 


and (9) , the necessary calculations being shown In columns 16, 17, 
an d 18. Actual values of h_ are obtained in column 19 by succes- 

11 hjv_;i 

slve multiplication by the ratio starting with the rim value 

% 


h^ determined from figure 4 for each value of (7 r ^ b and proceeding 

toward the center of the disk. The disk profiles and the stress 
distributions are shown in figure 5. The slenderness of the profiles 
relative to the rim is the result of operation at the elastic limit 
in most of the disk. More practical disks can be obtained by assum- 
ing o e values yielding greater factors of safety. The subject of 

allowable design stresses requires, however, considerable experi- 
mental investigation. 


SIMPLIFIED METHOD OF STRESS ANALYSIS OF GIVEN DISK 


Elastic Stresses 


A method for calculating the elastic-stress distribution in a 
disk of given design is presented in reference 3. Although this 
method can be directly used to check the stress distribution in any 
disk designed by the procedure described herein, simplifications can 
be made that appreciably reduce the amount of labor involved in 
calculation. 


Simplification of equations. - Brief ly, .the method of refer- 
ence 3 consists in rewriting stress equations (1) and (2) in finite 
difference form and solving these equations to obtain 


^r,n “ A r,n®t,a + ®r,n 

► 

(Ja * A. a. +B. 

VjQ W)I1 t;fl t,n J 


(14) 


where oi is the tangential stress at the first point station in 

the disk and remains unknown until the end of the calculation. The 
coefficients Ay n , and so forth, are determined from the equations 
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^n = K n A r,n-l +L n A t,n-l 
^,n 85 ^'n^n-l^'n^n-l 


B r,n = B n B r,n-l +I b B t,n-l' ,M n 


B* ts K' B ,+L' B^ ,+M' 
t,n n r,n-l n t,n-l n 






J 


(15) 


Thus, if the coefficients k T) A^, Bp, and B t , are known at 

the (n-l)st station, they are readily determined at the nth sta- 
tion from equations (15). The coefficients at the first station 
(r = a) are known. For a solid disk 


^,a = ^,a = 1 
B r,a 88 B t,a “ 0 

and for a disk with a central hole 

A r,a “ B r,a = ®t,a “ 

^t,a = 1 



(16) 


(17) 


Hence, successive applications of equations (15) yield the coef- 
ficients for subsequent substitution In equations (14) for actual 
stress determinations. The value of a Is determined from the 

condition at the rim where the radial stress Is known. Thus, from 
the first of equations (14) 


a 


t,a 


a 


r,b" B r,b 

Vb • 


( 18 ) 


where 0 r b is the radial rim stress and ^ and B r are the 

coefficients as determined In the last of the successive applica- 
tions of equations (15) . 

The bulk of the computational work in a stress analysis Is 
contained in the determination of the coefficients Kq, 1^, and 

so forth, in equations (15). When the stations are randomly chosen, 
the equations for these coefficients are complicated. A 



NACA IN 1957 


13 


great reduction In the amount of labor is 

are chosen according to an assigned value 

if Poisson’s ratio is assumed constant as 
slons assume the simplified form given in 


effected if the stations 

of n — . in addition, 
r n 

U = 0.333, the eipres- 
the following table: 


TABLE A 



r n r n-l = o.050 
*n 

r -r . 

n n-1 = 0.200 

Ee 

V 

0.959 - 0.007 

V ^n-l 

0.834 - 0.015 

V *n-l 


0.025 + 0.024 

V En-1 

0.104 + 0.077 

V Vl 

K ’n 

0.340 - 0.291 rr^- 

h n Eq-1 

0.343 - 0.153 

V V-l 

L *n 

0.009 + 0.942 

Vl 

lln — 1 

0.043 + 0.767 

V E n-1 

“n 

- ( 0 . 023^ip n _ j+0 . 025P n ) r Q 2 » 2 
-0.024E Q H’ n 

- (0 . 067^ip n _ 1+ O . 104P n )r fl 2 « 2 
-0.092E n H' n 

M ’n 

- (0 . 008^=ip n . 1+ 0 . 009p n )r n V 
-0.977E n H' n 

- ( 0 . 027^ip n _ 1 +0 . 043p n )r n 2 G) 2 

-0.920E H' 
n n 


Allnement charts. - Evaluation of the coefficients K^, K' n , 

and so forth from the expressions given in table A is possible. A 
further reduction in labor, however, can be effected by the use of 
the allnement charts in figure 6. For example, in figure 6(a), 

IL, 

joining with a straight line the points - = 0.987 and 

E n-1 

s 1.185 yields intersections on the suitable scales for values 
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Ljj = 0.199, L’ n = 0.808, ^ = 0.972, and K' n = 0.256, which 

check closely with calculations based on the expressions in table A. 
If the density can be considered constant between two stations 
(p n = P n _i)> the values of Mq and M' n can also be obtained from 

the charts. As an example, if fo-l „ 1.185, r n 2 po? « 20.5X10 3 , 

h 

and EqH'jj = 5.5X10, the corresponding values of — and 

r n W are joined by a straight line to determine the supplementary 
point P, which is then Joined to = 5.5x10^ to yield inter- 

sections Mjj = -4.3X10 3 and M' n = -6.6xl0 3 . Between any two sta- 
tions where a change in density occurs, the expressions for ^ and 
M' n in the table must be used. 

Illustrative example. - The necessary calculations are shown in 
table II for the determination of the elastic stresses in a disk of 
given design, which, are the disks formerly designed for the various 
assumed radial rim stresses. The coefficients K^, 1^, and so 

forth, listed in columns 20 to 25 are determined from values listed 
in table I used with the allnement charts in figure 6. Values beyond 
the range of the charts given in figure 6 may be computed from the 
equations given in table A. Columns 26 to 29 list the A and B 
coefficients as determined from equations (15) in which the values 
listed in equations (16) are used as the initial conditions. The 
value of C-fc a is then determined with the aid of equation (18), 

and the radial and tangential stresses at all stations are obtained 
from equations (14) and listed In columns 30 and 31, respectively. 

The octahedral sheen* stresses at various stations are listed in 
column 32, and the yield stresses at each station, which were the 
design criterions, are listed in column 33. A correlation of sat- 
isfactory engineering accuracy exists between the computed stresses 
and the design stresses shown in columns 14 and 15 of table I. 


Plastic Stresses 

Modification of elastic-stress equations. - A method for cal- 
culating the. stress distribution in a disk subjected to plastic 
flow is presented in reference 4. The form of the calculation is 
similar to that of an elastic-stress determination, except that the 
term H' n in the expression for M n and M' n (table A) is 

replaced by H' n -P ' n , where P' n depends upon the plastic strains 
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in the radial and tangential directions. In practice, the plastic 
strains are estimated and a stress distribution is determined. 
Comparison of the plastic strains based on the calculated stress 
values with the plastic strains assumed in order to determine these 
stresses permits better subsequent estimates. When resultant strain 
values agree closely with assumed strain values, the strains are 
assumed to be correct. 

The values of P' n to be used in the application of the sim- 
plified method are given in the following table: 


TABLE B 

r n" r n-l 

*n 

P n 

0.200 

0.100 A^ n + 0.125 A^n.! -1.100 A^ n + 0.875 A-^n-i 

0.050 

0.025 A r ^ n + 0.026 A r ^ n _^ -1.025 A t ^ n + 0.974 A^ n _^ 


where A_ is the plastic radial strain at the nth station, and 

the other A terms represent other plastic strains according to the 
notation of this report. The detailed procedure for estimating the 
plastic strains and checking with computed values is given in 
reference 4 . 

Similarly for creep calculations, Q' n terms are added to the 
P' n terms, the form of the expressions being identical to those 
for P ' in table B except that the A terms are replaced by creep 
terms & as discussed in reference 4. 

Illustrative example. - As an illustration of practical pro- 
cedure, the disk previously designed for a radial stress of 
30,000 pounds per square inch at the rim is considered. As pointed 
out in the design calculations, any redial stress above 0 pounds 
per square inch at the rim will cause the octahedral shear stress to 
exceed the elastic limit; hence, plastic flow must occur. The 
amount of plastic flow in this case is small and is confined to the 
single station at the rim, but the method involved is essentially 
the same for the more complicated cases involving plastic flow at 
several stations and previously accumulated plastic flow and creep. 
Examination of reference 4 will indicate the procedure in the more 
complicated cases . 
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The stress-strain curve for the disk material at the operating 
temperature of 1070° F is given in figure 7. Point A is the cal- 
culated elastic octahedral shear stress at the rim, B the yield 
point, and C the intersection with the stress-strain curve of a 
vertical line drawn through A. In accordance with the method of 
reference 4, the strain CD is the first approximation of the plastic 
strain. The value is listed in columns 34 and 35 of table III for 
station b. Columns 36 and 37 list the approximated radial and tan- 
gential plastic strains obtained from equations (20) of reference 4. 
The value of P' n , as determined from table B, is listed in col- 
umn 38. Column 39 corresponds to column 8 of table I used in the 
elastic-stress calculations; the only difference is that H' n -P ' n 

has replaced H' n . From this point, the calculations are identi- 
cal to elastic-stress calculations and are shown in columns 24, 25, 
and 28 to 32 (table III), which are similar to the corresponding 
columns of table II. Further approximations, made in accordance 
with the method of reference 4, are also shown in table III to 
obtain greater accuracy of approximation. 


COMMENT ON SIGNIFICANCE OF THEEMAL STRESSES 

The design method described assigns equal importance to the 
thermal and centrifugal stresses. The true relative importance of 
these two types of stress is, however, as yet to be evaluated. 
Whereas the centrifugal stresses are inherently necessary to bal- 
ance the centrifugal forces of rotation - so that yielding due to 
overstressing in one region shifts the load and thereby Increases 
the stress in another region - the thermal stresses constitute an 
internally balanced system. Tension in one region balances com- 
pression in others; hence, plastic flow in one region reduces the 
thermal stresses in another region. When considering bursting 
strength, for example, centrifugal stresses are probably of greater 
importance than thermal stresses. In the case of rim cracking 
(reference 8), which is in part due to cyclic plastic, flow at the 
rim, the thermal stresses are, however, fully as important as the 
centrifugal stresses. Likewise, when disk deformations are con- 
sidered in designing for close mechanical tolerances, thermal 
stresses must receive attention equal to that given to centrifugal 
stresses. 

The relative importance of centrifugal and thermal stresses 
depends, to a great extent, on the ductility of the disk material. 
As yet, experimental data are very meager on the significance of 
ductility. In any case, however, because the distinction between 
the two types of stress can only be realized after the occurrence 
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of plastic flow and because the behavior of a material In subsequent 
loading depends upon Its previous stress-strain history, thermal 
stresses should receive detailed attention early in the design In 
order to minimize any detrimental effects that may result from their 
omission. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, April 4, 1949. 
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APPENDIX - SYMBOLS 

In addition to the symbols used herein, the symbols of refer- 
ences 3 and 4 of which application is made in this report are defined. 
The various coefficients are given in their general form for random 
spacing of the point stations. From these coefficients, the sim- 
plified coefficients can be obtained by direct substitution of 

— — = constant and u = constant. 
r n 

A^ base area of airfoil section of blade, sq in. 

E elastic modulus of disk material, lb/sq in. 
h axial thickness of disk, in. 

N number of blades supported by disk 

r radial distance, in. 

T temperature, °F 

a coefficient of thermal expansion between actual temperature 
and temperature at zero thermal stress, in. /(in.) (°F) 

A plastic Increment of strain, in. /in. 

A r plastic Increment of strain in radial direction, in. /in. 

A.J. ' plastic Increment of strain in tangential direction, in. /in. 

AT temperature Increment above temperature of zero thermal 
stress, °F 

5 r creep increment in radial direction, in. /in. 

8^ creep increment in tangential direction, in. /in. 

C plastic strain corresponding to stress O q in tensile 
specimen, in. /in. 

p Poisson’s ratio 

p mass density of disk material, (lb) (sec^)/in.^ 
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0 stress at base of airfoil section of blade, lb/sq in. 

O octahedral shear stress, lb/sq in. 

O r radial stress, lb/sq in. 

a t tangential stress, lb/sq in. 

O v proportional elastic limit, lb/sq in. 

co angular velocity, radians/sec 

The following supplementary subscripts are used for denoting 
values of the preceding symbols in connection with the finite- 
difference solution: 


n nth point station 
n-1 (n-l)st point station 

a station at smallest disk radius considered (For disk with a 
central hole, this station is taken at the radius of the 
central hole; for a solid disk, this station is taken at a 
radius approximately 10 percent of the rim radius, which is 
larger than that (5 percent) recommended in reference 3. 

The value of 10 percent yields results of satisfactory 
accuracy and permits the use of fewer stations when a con- 
v -r i 

stant value of -£ — is used.) 
r n 

b station at rim of continuous section of disk or base of blades 

d station at rim of disk or base of airfoil section of blade 


The following supplementary symbols denote combinations of the 
previous symbols: 


V 


A 


stress coefficients defined by equations 




B. 


°r,n 


r,n 


B. 


l t,n J 


J t,n 


= A r,n°t,a +B r,n 
= A t,n CT t,a +B t,n 
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C n c * , a h n 


C 1 

Pn u. 

(l+UnKrn-rn-].) 

° n 

“ill 

2 E n r n 

Dn 

= |( r n 

■ r n-l^ 

D ’n 

n 

+ 

( 1+ ha> ( r n“ r n-l) 
2 E n r n 

F n 

- r n-lVl 

F ’n 

ha-1 

(l"H i n - 1 )( r n“ r n- 1 ) 

11 

1 

2 E n-l r n-l 

Gfc 

= I(r . 
2 V n 

■ r n-l)Vl 

p * 

. 1 

(l ^-l>< r xT r n-l> 

G n 

= Vi 

2 F n-l r n-l 


B n. = |" (r n- r n-l )(p nVn 2+p n-l h n-l r n-l 2 ) 

H’ n - VWA-l 

, _ r ’nWn 

^ C'nVCnD'n 

C n F 'n- c, n F n 

n 'O'A-V'a 

^nW’n 

^ " C'nW’n 
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L’ 


n 


_ CWn 

C'nWn 


Mn 


. P'nV^^'D-P'n^'n) 


C 'nW'n 


M' 


C'nVCntH'n-f'n-Q'n) 

C 'nWn 


(Mjj and M' n are defined in the elastic case for P' n = Q' n = 0) 


= A. 


r n:.^i 


r,n V 2r_ 


+ A. 


r,n-l 


- a* . i + Sa: g n-i 


2r, 


n-1 


t,n 


2r 


+ A* 


t > n-1 V ' 2 r n -i 


r n~ r n-l' 


Q* 


n 


6 r,nW7 + "*•“ 


Ka . », 


2r 


c f ! r n~ r n-l 1 

t>n ' 1 V 2r n .J 


For 


r n“ r n-l 


• = 0.050, 

x n 


E 


Z n = 1^1 (0 * 965 °t,n.r 0 - 298 °r,n-l)- E n H 'n 


r n _r n-l 

For - = 0.200, 

y% 9 


En 


Z n = ^ ^ (0.833 °t,n-l“ 9 * ^^n-l^-^n^'n 
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Thickness at base of serrated Temperature 

region, h h , in. - 
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Radial distance, r, in. 

figure 3. - Variation of temperature with disk 
radius assumed for illustrative design. 



Radial stress, a r 


lb/sq in. 


Figure 4. - Relation between thickness at 
rim of continuous section of disk and 
radial stress at that location. 
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